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In this paper, we carry out a theoretical analysis of the zero- frequency and finite-frequency shot 
noise in electron tunneling through a two-level interacting system connected to two leads, when a 
coherent coupling between the two levels is present, by means of recently developed bias-voltage 
and temperature dependent quantum rate equations. For this purpose, we generalize the traditional 
generation-recombination approach for shot noise of two-terminal tunneling devices properly to 
take into account the coherent superposition of different electronic states (quantum effects). As 
applications, analytical and numerical investigations have been given in detail for two cases: (1) 
electron tunneling through a quantum dot connected to ferromagnetic leads with intradot spin-flip 
scattering, and (2) spinless fermions tunneling through seriesly coupled quantum dots, focusing on 
the shot noise as functions of bias- voltage and frequency. 

PACS numbers: 72.70.+m, 73.23.Hk, 73.63.-b 
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I. INTRODUCTION 

In recent years the study of shot noise in mesoscopic 
quantum systems has become an emerging topic in meso- 
scopic physics, because measurement of shot noise can 
reveal more information of transport properties which 
are not available through the conductance measurement 
alonei^ For example, dynamic correlations in the tun- 
neling current originated from the Pauli exclusion princi- 
ple can provide information regarding the barrier geom- 
etry. On the other hand, strong Coulomb blockade effect 
also acts in correlating wavepackets and shot noise, when 
the charging energy becomes larger than the thermal en- 
ergy in small quantum devices. 

Besides, special interest has recently been put upon in- 
vestigation of tunneling through two- level system (TLS), 
when the coherent coupling between the two levels is pre- 
sented. In such internally coupled system, quantum in- 
terference effects between the two levels play a crucial 
role in its transport, inducing, for instance, temporal 
damped oscillations in average current and a correspond- 
ing fine structure in shot noise power spectrum S{lu) due 
to its intrinsic Rabi oscillationi2i^i^iSiL2i2ii£ Furthermore, 
it is predicted that a dip or a peak shown in S{uj) at 
the Rabi frequency reflects the differing relative phase of 
states carrying the tunneling current^ 

In literature, the quantum rate equations have been 
developed to describe this kind of quantum coherence 
effects in quantum tunneling through TLS, in associ- 
ated with the strong Coulomb blockade effecti ^^'^'^ So far, 
most calculations for shot noise in TLS have been based 
on the number-resolved version of quantum rate equa- 
tions, which takes numbers of tunneled electrons through 
TLS as parameters of density matrix elements (the degree 
of freedom of reservoirs), and have made use of Laplace 
transform to calculate the current-current correlation 
functions. Unfortunately, this scheme is only valid at 



zero temperature and high bias-voltage condition^SiLSili 
Therefore, it is quite desirable to find a way of evaluat- 
ing current correlation in arbitrary bias-voltages, partic- 
ularly in moderately small bias- voltage region, where the 
coherence plays a more prominent role in quantum tun- 
neling processes. Accordingly, the purpose of this paper 
is to develop a general scheme to study quantum shot 
noise spectrum of coupled small quantum systems based 
on the recently established quantum rate equations with 
a bias-voltage and temperature dependent version ►■''i 

Actually, a numerical method, named generation- 
recombination approach, was established for analysis of 
shot noise in two-terminal single-electron tunneling de- 
vices with the classical rate equation (classical shot noise) 
in Refs. (15lil6iil7i |. These earlier papers proposed a pow- 
erful method to evaluate the double-time current correla- 
tion function by switching the time-dependent current to 
a time evolution propagator of the density matrices in a 
consistent way. In the present paper we will modify this 
traditional approach to suit the underlying quantum rate 
equations for properly accounting for the quantum co- 
herence effects, i.e., the nondiagonal elements of density 
matrix, in calculation of quantum shot noise and develop 
a tractable computation technique in matrix form. 

The rest of this paper is organized as follows. In sec- 
tion II, we study the quantum shot noise for an interact- 
ing quantum dot (QD) connected to ferromagnetic leads 
with intradot spin-flip scattering. First we review the 
model Hamiltonian for this internally coupled TLS and 
the underlying quantum rate equations for description of 
coherent tunneling in sequential picture. Subsequently, 
we elaborate the general formalism for the quantum shot 
noise in matrix form and present a convenient numeri- 
cal method for the matrix calculation in frequency space 
with the help of the spectral decomposition technique. 
In the final part of this section, we analyze the zero- 
frequency Fano factor, which quantifies correlations of 
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shot noise with respect to its uncorrelated Poissonian 
value, in several limits, and then provide numerical re- 
sults for the Fano factor as functions of bias- voltage and 
frequency. In section III, we employ the same procedure 
to evaluate the quantum shot noise in two coupled quan- 
tum dots in series (CQD). Finally, a brief summary is 
presented in section IV. 

II. SHOT NOISE IN A SINGLE QUANTUM DOT 
WITH INTRADOT SPIN-FLIP SCATTERING 

A. Model and quantum rate equations 

The system that we study is a single QD with an ar- 
bitrary intradot Coulomb interaction U connected with 
two ferromagnetic leads. In this paper, we assume 
that the tunneling coupling between the dot and the 
leads is weak enough to guarantee no Kondo effect and 
that the QD is in the Coulomb blockade regime. This 
kind of spin-related single-electron devices suffers in- 
evitably intrinsic relaxations (decoherence) due to the 
spin-orbital interactioni or the hyperfine-mediated spin- 
flip transitions^ For simplicity, we model the intrinsic 
spin relaxation with a phenomenological spin-flip term 
R and assume that this spin-flip process happens solely 
within the dot (the spin is conserved during tunneling). 

Moreover, we assume that the temperature is low 
enough for one to see the effects due to the discrete 
charging and discrete structure of the energy levels, i.e. 
ksT ^ J7, A (A is the energy spacing between orbital 
levels). Each of the two leads is separately in thermal 
equilibrium with the chemical potential /i^, which is as- 
sumed to be zero in the equilibrium condition and is 
taken as the energy reference throughout the paper. In 
the nonequilibrium case, the chemical potentials of the 
leads differ by the applied bias-voltage V. We are inter- 
ested in the region eF <C A, where only one dot level 
{cd) contribute to the transport. Here we neglect Zee- 
man splitting of the level due to weak magnetic fields B 
(g^sB < ksT), which means that both the spin-up and 
spin-down currents through the dot go through the same 
orbital level e^. Therefore, the Hamiltonian of resonant 
tunneling through a single QD can be written as:^^ 

+4iCdT) + Vud^mx + Yl (V4fe<TCda + H-c), 

(1) 

where cj^^.^ {crfka) and c|j^ icda) are the creation (annihi- 
lation) operators for electrons with momentum fc, spin-a 
and energy trika in the lead (77 = L, R) and for a spin-a 

electron on the QD, respectively. Uda — c^da^da '^^ '^'^ 
cupation operator in the QD. The fourth term describes 
the Coulomb interaction among electrons on the QD. The 
fifth term represents the tunneling coupling between the 



QD and the reservoirs. We assume that the coupling 
strength V^^ is spin-dependent to take into account the 
ferromagnetic leads. 

Under the assumption of weak coupling between the 
QD and the leads, and with the application of the wide 
band limit in the two leads, electronic transport through 
this system in sequential regime can be described by the 
bias-voltage and temperature dependent quantum rate 
equations for the dynamical evolution of the density ma- 
trix elements*^ 

Poo = Yl^^'yP'^'^ ~ r^-rPoo), (2a) 

r\<j 

n n V 

+iR(pas - Psa), (2b) 



p„s = iR{Paa - pss) - ^ Y^^n-y + + + ^ris)P'ys, 

(2c) 

Pdd = E f ^iPTT +E r^^T/'ii -E(f .^T +^^i)Pdd, (2d) 

V n ri 

{a =t, i stands for electron spin and a is the spin oppo- 
site to <j). The statistical expectations of the diagonal 
elements of the density matrix, pu {i = {0, cr, d}), give 
the occupation probabilities of the resonant level in the 
QD as follows: poo denotes the occupation probability 
that central region is empty, p^a^ means that the QD 
is singly occupied by a spin-cr electron, and pdd stands 
for the double occupation by two electrons with differ- 
ent spins. Note that they must satisfy the normaliza- 
tion relation poo + pdd + J2a P<^<^ ~ ^- non-diagonal 
elements p^s describe the coherent superposition of dif- 
ferent spin states. These temperature-dependent tunnel- 
ing rates are defined as F^^,. = F^cr/^(ed) and F^^,. = 
Triaf^i^d + U)^ where F^o- are the tunneling constants, 
/+(cj) = {1 + e("~^")/^}^i is the Fermi distribution 
function of the 77 lead and /^(w) = 1 — f^{uj). Here, 
F^g. (F~g.) describes the tunneling rate of electrons with 
spin-(7 into (out of) the QD from (into) the 77 lead with- 
out the occupation of the a state. Similarly, F+^ (r^o-) 
describes the tunneling rate of electrons with spin-cr into 
(out of) the QD, when the QD is already occupied by an 
electron with spin-fj, exhibiting the modification of the 
corresponding rates due to the Coulomb repulsion. The 
particle current flowing from the lead 77 to the QD is 

^v/ e = Y'^^n<yPdd + ^^aP-yy - ^^sPaa - ^^aPoo)- (3) 

(7 

This formula demonstrates that the current is primar- 
ily determined by the diagonal elements of the density 
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matrix of the central region. However, the nondiagonal 
element of the density matrix is coupled with the diag- 
onal elements in the rate equation l)2b|l. and therefore 
indirectly influences the tunneling current. 



B. Quantum shot noise formula 

There is a well-established procedure, namely, the 
generation-recombination approach for multielectron 
channels, for the calculation of the noise power spec- 
trum based on the classical rate equations (classical shot 
noise) ^SiiSiil In this section we modify this approach in 
order to take into account the nondiagonal density matrix 
elements and derive the general expression for a quantum 
shot noise for the single QD. 

It is well known that the noise power spectra can be 
expressed as the Fourier transform of the current-current 
correlation function: 



Si^i^,{u:) = 2 / dte^-'[(/,(t)V(0))-(/,>(V)] 

^ — oo 

= 2(/,(t)V(0)).-2(/,)^(V)^. (4) 

Here, and I^i are the electrical currents across the 
and rj' junctions and t is the time. 

Before proceeding with investigation of current corre- 
lation, it is helpful to rewrite the quantum rate equations 
as matrix form: 



dp{t) 
dt 



M.p[t) 



(5) 



where p{t) = {pm. Pu^ Pii^ Pdd, Pm, Pi^Y is a vector 
whose components are the density matrix elements, and 
the 6x6 matrix M can be easily obtained from Eqs. ||2J). 
Therefore, the statistical averaging of an any time- 
dependent operator A{t) should be replaced by the sum- 
mation over all elements in this new p{t) representation: 



{A{t))^T,{Ap}^Y.^Apm 



(6) 



in which A is the matrix expression of the operator 
A under the basis {pm, Pn, Pa, Pdd, P^i, Pi\)'^ , and the 
summation goes over all vector [Ap(i)] elements (fc — 
1, 2, • • • ,6). Correspondingly, we can write the average 
electrical currents across the left [II) and right {Ib) junc- 
tions at time t as: 



{lL(R){t))^eY,^LiR)P{t) 



(7) 



where and Tji are current operators and the 

summation goes over all vector ^ L(R)P{t)] elements 
{k = 1,2,--- ,6). The current operators contain the 
rates for tunneling across the left and right junctions 
respectivelyfi^ and they can be read from Eq. Q as fol- 
lows: r~g. describes the tunneling rate of an electron with 



spin-cr out of the QD, when the QD were previously oc- 
cupied by two electrons. After this tunneling event, QD 
become singly occupied by a spin-a electron. In the other 
words, r~g. is a transition rate from pdd to the pss- It in- 
creases occupation probability p^g [this can be also seen 
by looking the second term of Eq. (|2b|) ] and decreases 
occupation probability pdd [the last term in Eq. (j2d|) ]. 
The term T^^p„ci in Eq. Q represents the tunneling of 
a electron out of the QD which leaves QD empty and 
can be interpreted as a transition rate from p„„ to 

the Poo ■ In similar way, T^g describes the transition rate 
from pgg to the pdd and is a transition rate from poo 

to the p„a- In the basis {pm, Pu, Pa, Pdd, P^i, Pi]V , the 
current operators have a matrix form: 



F^ = ± 



/ 





SI 







0\ 




-F+ 















-F+ 








St 





































\ 














0/ 



(8) 



where the sign of the current is chosen to be positive 
when the direction of the current is from left to right, so 
that the -I- sign in the last equation is for = L and the 
— sign stands for i] = R. The stationary current can be 
obtained as: 



TlP 



(0) 



Tb.p' 



(0) 



where is the steady state solution of Eq. lO 
which can be obtained from: 



(9) 



and 



Mp(") = 0, 



(10) 



along with the normalization relation ^^^j^ pl^-* = 1. We 
would like to point out that in our quantum version of 
rate equations, it is easy to check = (to = 

1,2,3,4), which implies that: 1) the Matrix M has a 
zero eigenvalue; 2) there is always a steady state solution 

3) the normalization relation X]n=i/'»(^) = 1 is 
independent on time. 

A convenient way to evaluate the double-time correla- 
tion function Eq is to define the propagator T{t) = 
exp[Mi], which governs the time evolution of the density 
matrix elements Pk(,t). The average value of the electrical 
currents across the left (II) and the right (Ir) junctions 
at a time t is given by: 



(/L(fl)(t)>=-e^[S(«)f(t)p(°) 



(11) 



which allows us to switch the time evolution from the 
vector p(t) to the current operators. Thus, we identify 
^L(R){'t) = ^L(R)T{t) as the time-dependent current op- 
erators. With these time-dependent operators we can 
calculate correlation functions of two current operators 
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taken at different moments in time. In particular, corre- 
lation function of the currents and in the tunnel 
junctions rj and 77', measured at the two times t and 
respectively, is given hy\i^ 



(/,(<)V(0)) = 0(<)^[f,(t)f„,p 



(0)1 



k 

6{-t)Y,[t,'f{~t)t^p^^\, (12) 



where 6{t) is the Heaviside function and the two terms in 
Eq. |{T^ stand for t > and for t <Q. The Fourier trans- 
form of propagator f{±t) is f{±uj) = {^iujl - 

where / is an unit matrix. We can further simplify this 
expression by using the spectral decomposition of the ma- 
trix M: 



M 



A 



(13) 



in which A are eigenvalues of the matrix M, S is a matrix 
whose columns are eigenvectors of M, E*^""^ is a matrix 
that has 1 at nn place and all other elements are zeros, 
and Px is a projector operator associated with the eigen- 
value A, so that r(±w) is 



T{±L0) 



Px 



A 



(14) 



Inserting expression for propagator T into Eq. H12|l 
current-current correlation in the w-space becomes 



(/,(i)V(0)>.. = 



E 

X,k 

E 

A,fc 



^riP\^r^' 
'■ \i 

—lUJ — A 



,(0) 



rV^\i\^(o) 

iw — \ 



(15) 



Eventually, substituting Eq. IjlSfl into the noise definition 
Eq. Q), and noting that in the summation over eigenval- 
ues the zero-eigenvalue contribution is canceled exactly 
by the term (/^)(/^/), we can obtain the final expression 
for a noise power spectrum: 




Sch 



Tri'PxTriP' 



(0) 



A 



lUJ - 



A 



(16) 



where S'^'^^ is the frequency-independent Schottky noise 
originating from the self-correlation of a given tunneling 



event with itself, which the double-time correlation func- 
tion Eq. (|12|l can not contain. Due to the fact that the 
current has no explicit dependence on the nondiagonal 
elements of the density matrix, it can be simply written 



as 
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Si 



Sch 



E 

k 



,(0) 



(17) 



The shot noise power spectrum is given by: 



S{uj) = a^Si^i^{uj) ^ l3'^Sij-,ij-,{uj) 

+al3Si,i,{oj) + al3Si,i,{uj). (18) 

Here, the coefficients a and /3, a + (3 = 1, depend on 
barriers geometry. -'^^ For simplicity we take a = (3 = 1/2. 
The Fano factor, which measures a deviation from the 
uncorrelated Poissonian noise, is defined as: 



F{u;) 



5H 
2el '' 



(19) 



where 2el is the Poissonian noise. 



C. Discussions and numerical calculations 

In the following we consider two magnetic configura- 
tions: the parallel (P), when the majority of electrons 
in both leads point in the same direction, chosen to be 
the electron spin-up state, a =t; Eind the antiparallel 
(AP), in which the magnetization of the right electrode 
is reversed. The ferromagnetism of the leads can be ac- 
counted for by means of polarization-dependent coupling 
constants. Thus, we set for the P alignment 

Plt = TflT = (1 + p)r, Tli = = (1 - p)r, (20) 

while for the AP-configuration we choose 

Tlt = r^i - (1 + P)T, Tli = TflT = (1 - p)r. (21) 

Here, P denotes the tunneling coupling between the 
QD and the leads without any internal magnetization, 
whereas p {0 < p < 1) stands for the polarization 
strength of the leads. We work in the wide band limit, 
i.e., P is supposed to be a constant, and we use it as an 
energy unit in the rest of this section. The zero of en- 
ergy is chosen to be the Fermi level of the leads in the 
equilibrium condition (/j.^ — /i^ = 0). The bias- voltage, 
V, between the source and the drain is considered to be 
applied symmetrically, fj,L = —fJ^R = eV/2. The shift of 
the discrete level due to the external bias is neglected. 

As a reference case for our analysis we use the analytic 
result for the case of paramagnetic electrodes, p = 0. In 
the system with paramagnetic electrodes, both channels 
for electrons with the spin up (|) and down (|) are equiv- 
alent and the tunnehng rates are J^r] — J2rj — 
and X),, ~ Srjf'f — We assume that Coulomb 
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repulsion is large and discuss the two special cases sepa- 
rately. First, we consider that no doubly occupied state 
is available in the QD, i.e, pdd = 0. In this case Fermi 
levels ^ of the two leads is bellow U+ed, meaning r+ ~ 0, 
r~ ~ r^CT- Then the quantum rate equations Eq. ^ 
become: 



factor are 



h 



TrT 



Ll R. 



2rL 



2rf 



(27) 
(28) 



Poo 
PTT 
Pll 

Pn 

■{i) 



r + r Pll- 2r+poo, 



-T-pii + r+poo + 2i?p 



Ripn - Pit) - r p^f^ 



(22a) 
(22b) 
(22c) 
(22d) 
(22e) 



where pj.^'' and pj^ stand for real and imaginary part 
of non-diagonal density matrix elements and we assume 
that the spin-flip term, R, is real. The matrix M is given 
as: 



M : 



F- 


r- 


-2T+ 








-F- 





F+ 





-2R 





-F- 


F+ 





2R 











-F- 





R 


-R 








-F- 



(23) 



The eigenvalues of the matrix M arc: Ao = 0,Ai = — F , 
A2 = -(F- + 2F+), X3 = XI = -{T- + 2iR) and the 
corresponding eigenvectors are given as the columns of 
matrix S: 



[(r-)2+2(r+)2]i/2 
T± 

[(r-)2+2(r+)2]i/2 


[(r-)2+2(r+)2]i/2 









1 


i 


i 


\ 


76 


V3 


^/3 





1 


i 


i 




76 


V3 


^/3 







2 

^/6 










1 



















1 


1 


/ 



(24) 



The current operators can be read from Eq. ^ and 
the nonzero elements are (rL)2i = — TJ, (rL)i2 = F^, 

(fL)3i = -rj, (fi)i3 - F^, and tn = -Fl- On 
the other hand, the steady states can be obtained from 
Eq. together with poo + J2a P<^<^ = 1, as: 



P?r = Pli 



F+ 



2F+ 



F- 



2F^ 



(25a) 
(25b) 



and the stationary tunneling current is II = Ir ~ 
P-^2r+ (-^Z^fl ~ r^F]^). Finally, applying the spectral 
decomposition of M and using Eq. (fTB|) one gets the Fano 
factor: 



F ^ 1 



4(rjF^ 



r^rj) 



2rjF£ 



(F- + 2F+)2 



^ R^ L 



^ R.^ L 



(26) 



In the case of large voltage, i.e., eT^/2 e^, F^ 
0, and FJ = F^, T~j^ = Tr. The current and the Fano 



F+ 

^ R 



The Fano factor depends only on the asymmetry in the 
coupling between the leads and the dot: it is equal to 
I for the completely symmetric case F^ = F^, and ap- 
proaches 1 when one of the coupling constants becomes 
much larger than the other one. 

In the opposite region, when the energies and ed + U 
are far below the Fermi level p, eV/2 Cd + U, we have 
F- = F+ = 0, fl = f+ ^ 0, and f+ ^ F+ = F,., 
F^ = F^ = Tji. Substituting these expressions into the 
rate equations and performing similar calculations, the 
electrical current and the Fano factor arc found to be 



I2 
F2 



FrF 



Li R 



= 1-2- 



FrF 



Li R 



(Fl 



(29) 
(30) 



The Fano factor is equal to ^ for completely symmetric 
coupling and to 1 for the asymmetric ones. 

When the leads are made of paramagnetic materials, 
the Fano factor does not depend on the spin-flip process 
and the same result was obtained by using the classical 
rate equationsii^ However, this is not true for the ferro- 
magnetic leads. 

Now we proceed our investigation of the Fano factor for 
the system where leads are ferromagnetic. Our numerical 
calculations for the current-voltage characteristic and the 
dependence of the Fano factor on bias- voltage for P- and 
AP-configurations are presented in Figs. 1-3. In these 
calculations we set = 1, Coulomb interaction U ~ -i, 
and temperature T = 0.1. By increasing the bias- voltage 
between two leads, two steps in the current-voltage char- 
acteristic occur: one is when the Fermi level of the source, 
PL, crosses the discrete levels Cd (for eV/2 > e^) and 
the other is when the Fermi level p^ crosses ed + U (for 
eV/2>€d + U). 

The effects of the polarization on the Fano factor with- 
out the spin- flip scattering {R = 0) are plotted in Fig. 1. 
An increase of the polarization will lead to an enhance- 
ment of the current noise in both configurations (P and 
AP) but for different reasons. Let us discuss the P- and 
the AP-configurations separately. 

When the leads are in the P-configuration [Fig. 1(a) 
and (c)] , an increase of the polarization will raise the tun- 
neling rates of electrons with the spin-up, F^i and F^i, 
but reduce the tunneling rates of spin-down electrons, 
Tli and Ffl|. Accordingly, this will induce an increase 
of the spin-up current and a decrease of the spin-down 
current but it will not affect the total current through the 
system as shown in Fig. 1(a), because the total current 
is equal to the summation of the spin-up and spin-down 
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FIG. 1: Current (a, b) and Fano factor (c, d) vs. bias-voltage 
in the P-configuration (a, c) and AP-configuration (b, d) cal- 
culated for different polarization p without spin-flip process. 
Other parameters are: td = 1, U = A, and T = 0.1. 



region, spin-up and spin-down electrons are tunneling 
through the different channels and there is no more com- 
petition between these two tunneling events. This leads 
to a reduction of the current fluctuation and the Fano 
factor becomes the same as in the paramagnetic case. 

The situation is completely different in the AP- 
configuration [Fig. 1(b) and (d)]. An increase of the 
polarization increases tunneling rates F^i and Fjj^ and 
decreases tunneling rates F i| and Tr-^ . An electron with 
the spin-up, which has tunneled from the left electrode 
into the QD, remains there for a long time because the 
tunneling rate F/j| is reduced by the polarization. This 
decreases the spin- up current. An increase of the polar- 
ization also decreases the spin-down current because it 
reduces the probability for tunneling of the spin-down 
electrons into the QD. This will decrease the total cur- 
rent through the system in the Coulomb blockade region, 
and at the same time increase the shot noise. For large 
voltage, in the region eV/2 > Cd + U, both conducting 
channels become available which results in reduction of 
the noise comparing with the Coulomb blockade region. 
In this case the Fano factor does not go to the paramag- 
netic value because the asymmetry in the tunneling rates 
are still presented. 



currents. In the case that the Coulomb interaction pre- 
vents a double occupancy of the dot, there will be com- 
petition between tunneling processes for electrons with 
up spin and those with down spin. The characteristic 
times for these two processes are unequal due to polariza- 
tion: there is fast tunneling for spin-up electrons but slow 
tunneling for spin-down electrons through the system. 
The electron spin which tunnel with a lower rate mod- 
ulates the tunneling of the other spin-direction electron 
(so-called dynamical spin-blockade).^'^ For a large value 
of polarization, it leads to an effective bunching of tun- 
neling events and, consequently, to the super-Poissonian 
shot noise. 

It is worth noting that our findings in the P- 
configuration seem to be in obvious conflict wiht the re- 
cent prediction in the sequential tunneling through multi- 
level systems, in which it was argued that a negative dif- 
ferential conductance (NDC) surely accompanied by the 
classical super-Poissonian noiseiSiiSS Actually, in these 
multi-level systems NDC can occur when two levels, sep- 
arated by nonzero energy difference Ae = £2 ^ ci, have 
different coupling strength to the leads. The tunnel- 
ing through the level which has a stronger coupling to 
the leads will be suppressed by the Coulomb interac- 
tion once the level with weaker coupling is occupied by 
electron. This reduces the total current, thus induces 
the NDC. However, this is not the situation in this pa- 
per. Here we assume small Zeeman splitting of the level 
(Ae = qhbB < ksT) and no energy difference between 
two spin levels. 

Further increasing the bias- voltage above the Coulomb 
blockade region, i.e, for eV/2 > Cd + U, opens one more 
conducting channel and removes spin-blockade. In this 




FIG. 2: Current (a, b) and Fano factor (c, d) in the P- 
configuration (a, c) and AP-configuration (b, d) calculated for 
polarization p — 0.5 and different spin-flip processes. Other 
parameters are the same as in Fig. Q 

The dependence of Fano factor on the spin-flip scat- 
tering with given polarizations p — 0.5 and 0.9 can 
be found from Figs. 2 and 3, respectively. Introduc- 
ing spin-flip scattering usually results in suppression of 
the zero-frequency current fluctuation and the Fano fac- 
tor. This behavior can be easily understood by the fol- 
lowing consideration. The spin-flip scattering opens ac- 
tually one new path for electrons to tunnel out of the 
QD: a spin-up(down) electron tunneling into the QD can 
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now experience spin-related interaction insider the dot 
and change its spin-polarization direction and then exit 
from the QD as a spin-down(up) electron. As a result, 
an electron which is forced to spend more time in QD 
due to polarized leads (for example, spin-down electron 
in the P-configuration and spin-up electron in the AP- 
configuration) now has greater probability of leaving the 
dot. In other word, the spin-flip scattering plays a com- 
pensating role in tunneling in contrast to that of the lead 
polarizations and consequently makes an opposite contri- 
bution to the current fluctuations. However, this is not 
true for the P-configuration when double occupation is 
allowed inside the QD (for eV jl > Ed + U). In this region 
the spin-flip scattering does not have any effect on the 
Fano factor [Figs. 2(c) and 3(c)]. Spin- up and spin-down 
electrons are passing through separate channels without 
changing their spins. 
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FIG. 4: Fano factor vs. polarization p in the P-configuration 
(a) and the AP-configuration (b) calculated for different spin- 
fiip processes R in the Coulomb blockade region, V — 4. 
Other parameters are the same as in Fig. Q 
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FIG. 3: The same as Fig. ^except for p = 0.9. 

Fig. 4 shows the Fano factor vs. polarization. In the 
P-configuration, the Fano factor increases with polariza- 
tion. In AP-configuration, for small spin-flip scattering, 
Fano factor increases with polarization but for larger R it 
starts to decrease. This is caused by competition of two 
effects: increase in the Fano factor due to polarization 
and, decrease due to spin-flip scattering. 

The frequency dependence of the Fano factor in the 
Coulomb blockade region is given in Fig. 5. For weaker 
spin-flip scattering, increasing scattering leads to reduc- 
tion in the Fano factor around zero-frequency, in both 
configurations [Fig. 5(a) and (c)]. While the larger values 
of the spin-flip scattering cause different low-frequency 
behaviors for different leads polarizations [Fig. 5(b) and 
(d)]. In addition, the strong spin- flip scattering also gen- 
erates an unambiguous peak in the P-configuration and 
a hump structure in the AP-configuration approximately 
located at a nonzero value of frequency u = 2R (the 



Rabi frequency), which characterizes the resonant oscil- 
lation between two spin states inside the QD. We provide 
a simple qualitative explanation for this behavior as fol- 
lowing: When spin-polarized electrons are injected into 
the QD from the left lead, the spin Rabi oscillation al- 
ways allows the electron to escape more easily from the 
dot to the right polarized lead, thus increasing the devi- 
ation of instant current from its average value, which 
induces enhancement or suppression of shot noise de- 
pending on the symmetry of states carrying currentpW In 
the P-configuration (Ft^i > F^jj^), the outgoing wavefunc- 
tion of electrons preserve the symmetry with the domi- 
nated incoming electrons in the left lead [because we set 
> for the polarization leads in this paper, see 
Eqs. (|20|l and H21() ]. which enhances the noise and thus 
causes a peak in the Fano factor [Fig. 5(a)]. Neverthe- 
less, the situation is complicated for the AP-configuration 
(Fj^i < Ffli). The anti-symmetry generates a hump 
structure and finally a dip in the shot noise spectrum 
for high spin-flip scattering rate [Fig. 5(c)]. Moreover, 
in the high frequency region, the Fano factor goes to a 
constant i for both configurations. 



III. SHOT NOISE IN COUPLED QUANTUM 
DOTS 



A. Model and quantum rate equations 

In this section we consider a resonant tunneling 
through a CQD with weak coupling between the QDs and 
the leads. We assume that electron hopping t between the 
two QDs is also weak. In order to simplify our derivation, 
we consider here the infinite intradot Coulomb repulsion 
[/' and a finite interdot Coulomb interaction U, which ex- 
cludes the state with two electrons in the same QD but 
two electrons can occupy different QDs. We are again in- 
terested in the regime kT <^ U, A, and eV A, so that 
only one level in each dot contribute to the transport. 
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FIG. 5: Fano factor in the P-configuration (a, b) and AP- 
configuration (c, d) vs. frequency calculated in the Coulomb 
blockade region for different spin-flip processes R. Other pa- 
rameters are: ~ 1, U = 4, V = 3, p = 0.5, and T = 0.1. 



-rjp22 - it{pi2 - P2i), (32c) 
Pi2 = +itipii - P22) 

-^[rz + r;^ + 5]f+]pi2, (32d) 

Pdd = FflPii + f +P22 - (r^ + f]^)pdd, (32e) 

with Poo + Pii + P22 + Pdd = 1- Here, poo denotes the 
occupation probabiUty that central region is empty, pu 
{i = 1,2) means that the ith QD is singly occupied by 
an electron, and pdd stands for the double occupation of 
the central region (each one of the QDs is occupied by 
one electron). The non-diagonal elements pij describe 
the superposition of the two levels in different QDs. The 
tunneling rates = r^/±(erf) and = r^/^(e<i + U) 
have the similar prescriptions as in the above section. 

The electric current II flowing from the lead L to the 
QD can be calculated as: 

Il/b = ripdd + r^pii - rjp22 - rjpoo- (33) 

Similarly, for the current flowing from the QD to the R 
lead we have 



The tunneling Hamiltonian for the CQD ii 



.14 



cr,cr 

E 

k.<7 



\0n2a' + ^^Lac\^,^C^^ + H.C.) 
(V'i?-4fe<xC2.+H.C.), (31) 



where cJ^2)ct' '^1(2)0- £^re creation and annihilation opera- 
tors for a spin-cr electron in the flrst (second) QD, respec- 
tively, ti (z = 1, 2) is the bare energy level of electrons in 
the ith QD. The other notations are the same as in the 
SQD case. 

Under the assumption of weak coupling between the 
QDs and the leads, and with the application of the wide 
band limit in the two leads, electronic transport through 
this system can be described by the master equation, 
Here, in order to simplify the analysis, we only consider 
spin independent tunneling processes and we take the 
bare mismatch between the two bare levels to be zero 
(ei = (-2 = Cd)- The desired quantum rate equations can 
be readily obtained as: 



Poo 



TiPii + ^rPi2 ~ (rj + r+)po 



pii = rjpoo + r^prfd - r^pii 
-r;^pii + i<(pi2 - P21), 

P22 = r^poo + r^prfd - r^p22 



(32a) 
(32b) 



^fl/e = -(r^pdrf + r^p22 - r^pii - Tj^poo). (34) 

Again, the sign of the current is chosen to be positive 
when the direction of the current is from the left to the 
right. It is easy to prove that, in stationary condition, 
the current conservation is fulflUed Ij^ = Ir- 



B. Discussions and calculations of shot noise 

In order to calculate the noise power spectrum in CQD, 
we can simply employ the same procedure that we de- 
scribed in Sec. II. First of all, the simplified quantum rate 
equations Eq. (|32|l can be rewritten in the matrix form 
Eq. (O, in which p(t) = (pao, Pu, P22, Pdd, P12, P21V and 
the matrix M can be read from Eq. (|32|l . Next, the cur- 
rent operators can be obtained from Eqs. H33|) and l|34|) 
as follows: (F^) describes a tunneling of an electron 
from the first (second) QD to the L(i?)-lead and it can 
be considered as a transition rate from pdd to P22 (pii)- 
Similarly, F^ {^r) is a transition rate from pn (P22) to 
Poo, (r^j) describes transition from P22 (pii) to pdd, 
and FJ (Fjj) is a transition rate from poo to pn (p22)- In 
the basis {poo, Pii, P22, Pdd, P12, P21V , the left and right 
current operators are 
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(35) 
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and 
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(36) 



respectively. Finally, we can make corresponding substi- 
tution in Eqs. i{TU\i . lfTC|l . and ((T7|l to compute the quan- 
tum shot noise in CQD and use Eq. p9|l to calculate the 
Fano factor. 

In the following discussions, we set parameters of the 
CQD under consideration as: = 1, Coulomb interac- 
tion U — 4, and temperature T ~ 0.1. And we also 
apply the bias voltage V between two leads symmetri- 
cally, eV/2 = fj,L = —fJ.R- The zero of energy is chosen 
to be the Fermi level of the leads in the equilibrium and 
the energy unit is the one of the tunneling constants. 

First, we analyze the zero- frequency shot noise in three 
limits: (i) zero voltage limit; (ii) the Coulomb block- 
ade region; and (iii) high voltage limit. When the bias 
voltage is below the resonance, i.e., eV/2 < e^, the 
transport through the system is not energetically al- 
lowed and the noise is Poissonian {F = 1). As the 
bias voltage is increasing, there occur consecutively two 
plateaus in the current-voltage characteristic, separated 
by the thermally broadened step, corresponding to the 
case of eV/2 > ej, (the Coulomb blockade region) and of 
eV/2 > Ed + U (high voltage condition), respectively. On 
the first plateau, where only one electron can be allowed 
inside the system, the current Ii and the Fano factor Fi 
become: 



h 



R^ L 



^'r 



8^2 



[4t^2^L + ^R) + TLrlf' 



(37) 



(38) 



While in the second plateau, the CQD is in doubly oc- 
cupied state, and the current I2 and the Fano factor F2 
are: 



{Tl + Tr) (TlTr + 4t2) ' 



F2 = l- STLTRt^ 



STlTr + TI + TI+ Ae 
[TL + TRf [TLTR + U^ f 



(39) 



(40) 



These results are in agreement with the previous deriva- 
tion from the Laplace transform in Ref. [g ■ 

We proceed with numerical calculations. In Fig. 6, we 
plot the currents and zero-frequency Fano factors in the 
CQDs as a function of bias voltage for various hopping 
rates between dots t/T = 0.2, 0.5, 0.7, and 1.0. We 



find that, depending on the relation between the hop- 
ping t and the coupling constants Yr{l) ■, the CQD man- 
ifests three distinct physical scenarios in its tunneling 
and low- frequency fluctuation properties: (i) In the case 
of hopping 2t > ^/TlTr {t > 1/2 in Fig. 7), current 
increases in the second step {I2 > h) and the system 
exhibits positive differential conductance (PDC). Corre- 
spondingly, the Fano factor decreases {F2 < Fi) with 
increasing bias-voltage; (ii) If the hopping 2t = y/TLTR 
{t — 1/2 here), two currents and shot noise are equal to 
each other, respectively, at the entire region; (iii) While 
for the hopping 2t < ^/^L^R {t < 1/2 here), current de- 
clines in the second step and consequently negative differ- 
ential conductance (NDC) appears, but the Fano factor 
raises. Moreover, we find that the Fano factor is always 
below 1 in all three regimes, indicating that the quantum 
shot noise is perpetually smaller than the classical value 
(sub-Poissonian). Once again, we obtain results differing 
from those by the recent papers Refs. ^21 22]: Our shot 
noise is enhanced but does not reach super-Poissonian 
value in the NDC regime. 
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FIG. 6: The calculated current (a) and Fano factor (b) vs. 
bias voltage in the CQDs with different hoppings t between 
two QDs. Other parameters are: Fl = r_R = F, = 1, 
[/ = 4, and r = 0.1. 

Fig. 7 and 8 shows the current and the Fano factor as 
a function of the tunneling constants F^ and F^ in the 
current plateaus regions. For small hopping t <Vl the 
currents I\ and I2 show non-monotonic behavior with in- 
creasing Y R [Fig. 7(a)]: They increase for small couplings 
and, decrease when couplings are larger. With increasing 
Fi [Fig. 7(b)], I\ increases for small couplings and grad- 
ually saturates. On the contrary, I2 shows an analogous 
non-monotonic behavior as in Fig. 7(a). The NDC ap- 
pears whenever I2 becomes smaller than 1\ and it is most 
pronounced for larger coupling to the left lead F^ ^ Yr. 
In small hopping limit (t < Y^^r)), Eqs. and (jSHJ 
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can be simplified, respectively, to: 
and 



At' 



r. 



(41) 



(42) 



where Td^ = ^i^ and F^^ = Fi^ + Fl are the decoherence 
rates due to the interaction with the leads. Actually, if we 
simplify Eq. (|32d|l in the regions of the first and second 
current plateaus, we obtain: 



Pi2 = it[pii - P22) - ^rdiPi2, 



and 



P12 = it{pii - P22) 



(43) 



(44) 



respectively. It is clear that the two rates indeed charac- 
terize the decay of coherent superposition. Conservation 
of current allows us to express the current through the 
system as / = it{pi2 — P2i)- From Eqs. (|43|l and (|44|l one 
can see that for large decoherence rates, F^^^j ^ t, the 
probability of tunneling between the two dots, pi2, be- 
comes purely damped with time, implying the destructive 
quantum interference between the coherent superposition 
of two quantum dot states and the electronic state in 
the continuum (right lead) in the process of tunneling^S^ 
which induces localization of the electron in the first dot 
(QDs form ionic- like bonds), as well as the reduction of 
the two resonant currents Ii and l2- 

In the case of high voltage, the presence of one excess 
electron inside the dots further destroys coherent super- 
position because there is no possibility of electron tun- 
neling between the two dots due to the infinite intradot 
Coulomb interaction. The decoherence rate in this case 
is a sum of the two rates: F^, which describes the tun- 
neling of the second electron inside the system and Tn, 
which stands for the tunneling of an electron to the right 
lead. The decrease of the current I2 is caused by both of 
these processes [Figs. 7(a) and (b)], therefore it is smaller 
that Ji and NDC appears. In the NDC regime the noise 
increases but it remains sub-Poissonian (Fig. 8). 

In the large hopping limit [t 3> rdi(2))) electron 
can shuttle between the two dots many times in a phase- 
coherent way before it tunnels out into the right lead 
and thus becomes delocalized (covalent-like bonds). In 
this limit, the current and the Fano factor approach to 
the single QD values Eqs. lf77|l - (|Sn|l . 

The Fano factor is plotted against normalized fre- 
quency in Fig. 9. Similar noise properties are found 
for a coherently coupled double well structurei^ In the 
Coulomb blockade regime [Fig. 9], if an electron from 
the left lead is injected into the first QD, no further elec- 
trons can enter in the QD until this electron is removed. 
The time scale for this removal process is determined by 
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FIG. 7: Current vs. bias-voltage in the CQDs calculated for 
different coupling constants r_R (a, c) and Fl (b, d). Other 
parameters are: U — 4, ea = 1, T — 0.1, t — 0.5 (a, b) and 
t^2 {c, d). 
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FIG. 8: Fano factor vs. bias-voltage in the CQDs calculated 
for different coupling constants Fl (a, c) and Vr (b, d). Other 
parameters are as in Fig. |7| 



t^^ . Thus, in small hopping limit, the zero- frequency 
shot noise is reduced with increasing t due to the inter- 
dot Coulomb blockade effect [Fig. 9(a)]. As the frequency 
increases, the electron inside the first QD has more op- 
portunity to instantly tunnel into the second QD, which 
enhances the shot noise. While for the high values of 
hopping, the electron in the second QD can either escape 
to the right lead, which takes place on a time scale deter- 
mined by F^^, or, it can periodically return to the first 
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QD at a frequency 2t (the Rabi frequency). If 2t ^ Tji, 
there is larger probability for the electron in the second 
QD to tunnel back into the first QD than to escape to the 
right lead, which prevent another electron from entering 
the first QD. Thus at large values of t noise suppression 
occurs at a; = 2t, the dip in the Fano factor [Figs. 9(b)]. 




FIG. 9: Fano factor vs. frequency calculated in the Coulomb 
blockade regime for different hoppings t between the QDs. 
Other parameters are: Tr = Tl=^,V = 6, ed = l,U = 4, 
and T = 0.1. 

The influences of the coupling constants on the 
frequency-dependent noise characteristic is given in 
Fig. 10. When only one electron is present inside the 
CQD system [Figs. 10(a) and (b)], decreasing coupling 
to the right lead (decreasing decoherence rate) leads to 
the noise enhancement in the low frequency region. On 
the other hand, the decreasing decoherence rate results 
in the transformation of the electronic states in the CQD 
from the ionic-like bond to the covalent-like bond and the 
formation of a dip in the shot noise spectrum at moder- 
ately high frequency co = 2t. In the large voltage region, 
when doubly occupancy is allowed, the decoherence rate 
is given as a sum ofTL and Tji and both coupling con- 
stants influence shot noise in the same way [Figs. 10(c) 
and (d)]. 

IV. CONCLUSION 

In this paper we have presented theoretical investi- 
gations of the shot noise spectrum in resonant tunnel- 
ing through a interacting quantum coupled system, in 
which quantum interference effects play an important 
role in its transport and fluctuation properties. For this 
purpose, we have modified the well-known generation- 
recombination approach, which has been established for 
over ten years to study the shot noise in single-electron 
tunneling devices based on the classical rate equations, 
to incorporate with the quantum version of rate equa- 
tions. We have also developed a convenient numerical 
technique to compute the shot noise by applying the ma- 
trix spectral decomposition. 

As applications of our formalism, we have systemati- 
cally analyzed the current shot noise, as functions of bias- 
voltage and frequency, through (i) a single QD connected 
to two ferromagnetic leads with weak intradot spin-flip 
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FIG. 10: Fano factor vs. frequency calculated in the Coulomb 
blockade regime (a, b) and in the largo voltage region (c, d) 
for different coupling constants Fl and Tr between the QDs 
and the leads. Other parameters are: t = 0.5, Tr = T (a, c), 
Tl = T {h, d) V = 6 (a, b), V = 12 (c, d), ea = 1, U = 4, 
and T = 0.1. 



scattering, and (ii) coherently coupled QDs. The influ- 
ence of the Coulomb interaction and coherent electron 
evolution on the Fano factor has been investigated in 
detail. First we have given some analytic expressions 
for the zero-frequency Fano factors in two special cases: 
Coulomb blockade regime and double occupation region 
(high bias- voltage limit). It is shown that these results 
are in perfect agreement with previous analysis in the 
literature. 

In addition, we have performed numerical simulations 
for the two systems. For the single QD, we found 
that: 1) In the Coulomb blockade regime, an increase 
of the polarization leads to an enhancement of the zero- 
frequency current noise in both configurations, even a 
super-Poissonian value in the P-configuration due to dy- 
namical spin-blockade, implying bunching effect in tun- 
neling processes. For the AP-configuration, this enhance- 
ment results from the asymmetry in the tunneling rates 
into and out of the QD (Fli > F^i but Tn < F^i) for 
each spin separately; 2) The spin-flip scattering compen- 
sates for polarization of the leads (removing the dynam- 
ical spin-blockade) and causes reduction of the current 
fluctuation and Fano factor; 3) The frequency-dependent 
Fano factors clearly show a peak in the P-configuration 
but a dip in the AP-conflguration at the Rabi frequency 
ui = 2R, reflecting differing symmetries of states carrying 
current. For the CQD, we have predicted that depending 
on the relation between the hopping t and the dot-leads 
couplings rj;,(/{), there are three distinct physical scenar- 
ios in its tunneling and low-frequency fluctuation prop- 
erties. More importantly, we found appearance of NDC 
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in the current-voltage characteristic and corresponding 
enhancement of the Fano factor but still remaining sub- 
Poissonian in the case of < < \/TiXb . Our sequent quali- 
tative analysis claimed that this feature is a result of the 
enhanced decoherence rate induced by the presence of the 
second electron inside the system [Eqs. (|41f) and H42|) ]. 
Moreover, the destructive interference between two QDs 
states leads to a dip in the frequency-related Fano fac- 
tor at the Rabi frequency 2t, which is also controlled by 
the tunnel coupling to the right lead in the Coulomb 
blockade regime. 

Expectably, our results show that shot noise measure- 
ments can provide more useful information about quan- 



tum coherence of electron wavefunction than do conven- 
tional transport measurements. 
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